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ABSTRACT 

We have carried out an investigation of the abundance of deuterium along two extended sight lines 
through the interstellar medium (ISM) of the Galactic disk. The data include Far Ultraviolet Spectro- 
scopic Explorer (FUSE) observations of HD 195965 (BILb) and HD 191877 (B0V), as well as Space 
Telescope Imaging Spectrograph (STIS) observations of HD 195965. The distances to HD 195965 
and HD 191877, derived from spectroscopic parallax, are 794 ± 200 pc and 2200 ± 550 pc, respec- 
tively, making these the longest Galactic disk sight lines in which deuterium has been investigated 
with FUSE. The FUSE spectra contain all of the H I Lyman series transitions (and the correspond- 
ing D transitions) except Lya. The higher Lyman lines clearly show the presence of deuterium. We 
use a combination of curve of growth analyses and line profile fitting to determine the D I abun- 
dance toward each object. We also present column densities for O I and N I toward both stars, 
and H I measured from Lya absorption in the STIS spectrum of HD 195965. Toward HD 195965 
we find D/H=(0.85±g'^) X 10" 5 (2a), 0/H=(6.61±}??) x 10~ 4 , and N/H=(7.94±^) x 10" 5 . Toward 
HD 191877 we find D/H=(0.78±^§) x 10~ 5 (2a) and N/H=(6.76±?;p) x 10~ 5 . The O I column density 
toward HD 191877 is very uncertain. Our preferred value gives O/H=(3.09±g g|) x 10~ 4 , but we cannot 
rule out O/H values as low as 0/H=1.86 x 10~ 4 , so the O/H value for this sight line should be taken with 
caution. The D/H ratios along these sight lines are lower than the average value of (1.52 ± 0.15) x 10~ 5 
(2a in the mean) found with FUSE for the local interstellar medium (~ 37 to 179 pc from the Sun). 
These observations lend support to earlier detections of variation in D /H over distances greater than a 
few hundred pc. The O/H ratio toward HD 195965 is supersolar. This star is part of an OB associa- 
tion, so there may be local enrichment by nearby massive stars. The D/H and O/H values measured 
along these sight lines support the expectation that the ISM is not well mixed on distances of ~ 1000 
pc. These observations demonstrate that although D/H studies through Lyman absorption may become 
impractical at d>2500 pc and log N(H I)>21, D/H studies in the distance range from 500 to 2500 pc 
may be very useful for investigating mixing and chemical evolution in the ISM. 

Subject headings: ISM: abundances — ultraviolet: ISM — stars: individual (HD 195965, HD 191877) 



1. INTRODUCTION 

The abundance of deuterium in the interstellar medium 
(ISM) is a strong constraint on chemical evolution sce- 
narios. Unlike most elements, the net cosmic abundance 
of deuterium is not enhanced by stellar processes or su- 
pernovae, so its abundance has been declining from the 
primordial value due to destruction by stellar processing. 
The primordial abundance of deuterium is itself an impor- 
tant quantity, as it can be used to measure r], the univer- 
sal ratio of baryons to photons, which can be compared 
to the predictions of Big Bang nucleosynthesis (Reeves ct 
al. 1973; Boesgaard & Steigman 1985; Buries, Nollett, & 
Turner 2001). 

One approach to finding the initial deuterium abun- 
dance is to measure it in high-rcdshift QSO absorption- 
line systems, assuming that the clouds which give rise to 
these systems have undergone very little stellar processing 
(Webb et al. 1997; Tytler et al. 1999; O'Meara et al. 2001; 
Pettini & Bowen 2001; Levshakov et al. 2002). However, 
the large range in D/H values measured in these systems 



suggests that the gas in some of these clouds may in fact 
have undergone some stellar processing (Prantzos & Ishi- 
maru 2001; Fields et al. 2001). It is therefore crucial to 
understand how chemical evolution affects the abundance 
of deuterium relative to other elements. Several groups are 
measuring the abundance of deuterium and other elements 
along many sight lines in the Galactic ISM in order to es- 
tablish the mean D/H ratio in the Galactic disk, and to 
look for variations in the D/H ratio. These variations can 
shed light on the effects of star-formation history and mix- 
ing processes on the D/H ratio. The local D/H value can 
then be compared to regions that have undergone varying 
amounts of stellar processing (e.g., high velocity clouds, or 
intergalactic H I clouds), which will illuminate the effects 
of stellar processing on the D/H ratio in these systems. 

The Copernicus satellite was used to measure the D I 
and H I Lyman series absorption toward ~ 20 stars, and 
the results suggested a D/H ratio of 1.5xl0 -5 , with a dis- 
persion that may reflect real variations in the ISM (see 
Vidal-Madjar & Gry 1984 for a review). Several measure- 
ments of D/H in the local ISM have been made with the 
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Hubble Space Telescope (HST) (Linsky et al. 1995; Linsky 
& Wood 1996; Linsky 2003; Lemoine et al. 1996; Vidal- 
Madjar et al. 1998; Hebrard et al. 1999; Sahu et al. 1999; 
Venncs et al. 2000), though only the Lya transitions of 
D I and H I could be analyzed, limiting the number of 
available sight lines to those with low enough H I columns 
that the D I absorption is not blended with the H I line 
but large enough H I columns that D I can be detected. 
Three sight lines observed with the Interstellar Medium 
Absorption Profile Spectrograph (IMAPS) spanned larger 
distances than the HST sight lines (250 to 500 pc versus 
< 100 pc), and like the Copernicus results they showed 
evidence of variation in D/H with location (Jenkins et al. 
1999; Sonneborn et al. 2000). 

One of the science goals of the Far Ultraviolet Spectro- 
scopic Explorer (FUSE) is to measure D/H in the Milky 
Way (Moos et al. 2000). FUSE is uniquely suited to do 
this because it observes the spectral region from 905 to 
1187 A which includes all of the D I Lyman transitions 
except Lya. The ability to use multiple D I transitions 
means that FUSE is able to measure 7V(D I) along higher 
column density sight lines than missions that only ob- 
serve Lya (e.g., HST), and for many more sight lines than 
short term missions (e.g., IMAPS). Moos et al. (2002, 
and references therein) reported results for the first set of 
seven sight lines ranging from 37 to 179 pc, five of which 
had both D I measurements (from FUSE) and H I (from 
HST). The results were consistent with a single value of 
D/H= (1.52 ± 0.08) x 10 -5 (the uncertainty is lcr in the 
mean). However, none of these sight lines extend as far as 
the three sight lines observed by IMAPS. Analysis of more 
extended sight lines is crucial for the verification and un- 
derstanding of variations in the D/H ratio in the ISM. 

In this paper we present FUSE measurements of D/H 
toward two distant stars in the Galactic disk: HD 195965 
and HD 191877. These stars are at distances greater than 
the IMAPS targets, extending the distance and range of 
N(H I) for D/H measurements in the ISM. Information 
about these two sight lines is listed in Table 1. The FUSE, 
STIS, and ground-based observations and data reduction 
are described in § 2. Section 3 describes some properties 
of the sight lines. In § 4 we describe both the curve of 
growth and the absorption profile fitting analyses of the 
D I column density toward both stars. In § 5 we discuss 
the determination of the O I and N I abundances along 
these sight lines, and in § 6 we re-examine the H I column 
densities. Section 7 contains a discussion of the results 
and conclusions. 

2. THE DATA 

2.1. FUSE Observations 

Information for the FUSE observations is given in Ta- 
ble 2. The FUSE mission and instrument are described by 
Moos et al. (2000) and Sahnow et al. (2000). The spectra 
were obtained through the large (30"x30") LWRS aper- 
ture. FUSE consists of 4 co-aligned optical channels, two 
optimized for longer wavelength (LiFl and LiF2; 1000 - 
1187 A) and two optimized for shorter wavelengths (SiCl 
and SiC2; 905 - 1100 A). All of the deuterium lines used 
in this paper are located in the SiC channels. The Ly/3 
D I line was not detected in either spectrum because it 
was blended with the strong H I line, but some lines of 



other elements used in our analysis are located in the LiF 
channels. Between exposures the focal plane assemblies 
were moved so that the spectra sampled different portions 
of the detector. This shifting lessens the effects of detector 
defects when the individual spectra are combined. 

The raw spectra were processed through the FUSE cali- 
bration pipeline (CALFUSE v2.1.6) at the Johns Hopkins 
University. The pipeline screens data for passage through 
the South Atlantic Anomaly and low Earth limb angle 
pointings and corrects for thermal drift of the gratings, 
thermally-induced changes in the detector read-out cir- 
cuitry, and Doppler shifts due to the orbital motion of the 
satellite. Finally, the pipeline subtracts a constant detec- 
tor background and applies wavelength and flux calibra- 
tion. After the pipeline reduction the individual spectra 
were co-added to produce the final calibrated spectrum. 
The data from the four channels were analyzed separately, 
because there are slight differences in the spectral resolu- 
tion between channels, which can cause problems if the 
channels are co-added. This has the benefit of allowing at 
least two measurements of each line, a further safeguard 
against detector defects. 

The presence of bright airglow lines at specific regions of 
the spectrum has resulted in a long-term depletion of elec- 
trons in the detector in those regions, and this can cause 
localized shifts in the spectrum in the dispersion direction. 
This effect, called "x-walk," may change the shapes and 
areas of absorption lines. The effect is more pronounced 
in certain FUSE detectors. We have avoided using the 
equivalent widths of absorption lines in regions where the 
x-walk effect is pronounced, so these detector artifacts do 
not affect the analysis in this paper. The x-walk distor- 
tion is described in the FUSE Observer's Guide available 
at the FUSE website (http://fuse.pha.jhu.edu/). 

Figure 1 shows a portion of the FUSE spectrum of 
HD 195965, and Figure 2 shows the same spectral region 
for HD 191877. At A - 920 A, the shortest wavelengths 
used in this paper, the signal to noise ratio is >30: 1 per 
~20 kms -1 resolution element in both spectra, and is 
considerably higher at longer wavelengths. 

2.2. STIS Observations 

We observed HD 195965 with the Space Telescope Imag- 
ing Spectrograph (STIS) on the HST on 2001 October 9 as 
part of observing program GTO-8487 (P.I.: Moos). Brief 
information about the observation is given in Table 2. 
The spectrum was obtained in ACCUM mode using the 
FUV MAMA, the E140H grating (A c = 1307 A), and the 
0.1" x 0.03" aperture. We reduced the data using CAL- 
STIS (v2.12a) with all of the normal reduction procedures 
enabled in the standard data pipeline, including the two- 
dimensional scattered light removal algorithms (see Lands- 
man & Bowers 1997; Lindler & Bowers 2000; see also Howk 
& Sembach 2000). A description of the STIS pipeline pro- 
cessing steps can be found in the HST Data Handbook 
(Brown et al. 2002). 

The extracted STIS spectrum has S/N - 25 - 30 per 
spectral resolution element. The spectral resolution of 
~ 2.6 kms" 1 (R = A/AA - 200,000) is « 11 times bet- 
ter than that available with FU SE. The absolute wave- 
length scale is accurate to 0.5 — 1.0 pixel, or about 3-6 
mA (~ 0.6 - 1.3 kms" 1 ; see Leitherer et al. 2001). The 
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post-processed STIS data have very low residual back- 
grounds. We find that the strongly saturated cores of the 
O I 1302.169 A and C II 1334.532 A lines have residual 
intensities of < 1%. 

2.3. Ground-based Observations 

High-resolution echellc spectra of HD 195965 and 
HD 191877 were taken with the 2.7 m telescope of the 
University of Texas McDonald Observatory on 2001 De- 
cember 29. The "2dcoude" spectrograph (Tull et al. 1995) 
was used with 1 camera setup and a Tektronix charge- 
coupled device (CCD) as the detector, which provided 
high spectral resolution (R = A/AA « 160,000; Av « 
2 km s" 1 ) and covered the wavelength range 6000 A to 
8600 A. This high resolution mode yielded spectra con- 
taining 12 disjoint echelle orders approximately 40 A wide, 
separated by approximately 150 A at 6000 A. The inter- 
order separation increased with wavelength until there was 
approximately 280 A between orders at near infrared wave- 
lengths. This setup provided spectra on interstellar K I. 
The resolution was determined from the full width at half 
maximum (FWHM) of Thorium- Argon (Th-Ar) calibra- 
tion lamp lines. 

The data were reduced in a standard way utilizing the 
NOAO SUN/IRAF software (v2.11.3). Dark, bias, and flat 
field lamp exposures were taken to remove any instrumen- 
tal effects due to the CCD detector. Comparison spectra 
were taken periodically throughout the night, typically ev- 
ery two hours. The average bias exposure was subtracted 
from all raw stellar, comparison (Th-Ar) and flat images. 
The scattered light was fitted by low order polynomials in 
both the dispersion and spatial directions, and removed. 
Pixel to pixel sensitivity variations were removed by divid- 
ing the normalized average flat field into the stellar spec- 
tra. The extracted spectra were placed on a heliocentric 
wavelength scale using the Th-Ar comparison spectra and 
corrections for the Earth's motion. The spectra were co- 
added and normalized to unity yielding a final spectrum 
with high signal to noise, approximately 100:1. Figure 3 
shows the spectral region near the K I 7698.964 A line for 
both sight lines after correcting the velocities to the Local 
Standard of Rest (for HD 195965 V LSR = V heUo + 16.7 
kms- 1 , and for HD 191877 V LSR = V heHo + 17 .4 kms" 1 ). 

3. PROPERTIES OF THE SIGHT LINES 

HD 195965 (Z=85.71°, 6=5.00°) was classified by Mor- 
gan, Code, & Whitford (1955) as type B0V. It is listed 
by Humphreys (1978) as belonging to the Cygnus OB7 as- 
sociation, which is at a distance of ~830 pc in the Local 
spiral arm. Diplas & Savage (1994) give a distance of 
794 pc from spectroscopic parallax (with an uncertainty 
of ~ 25%). At this distance it would be 69 pc above the 
Galactic plane, still well within the thin disk, where most 
of the cool gas is located. The Hipparcos catalog lists 
an uncertain distance of the star at d= 520+^0 pc (ESA 
1997). Diplas & Savage (1994) find log iV(HI)=20.9 based 
on IUE data. 

HD 191877 (Z=61.57°, 6=-6.45°) is type Bl.OIb (Mor- 
gan et al. 1955), and spectroscopic parallax gives a dis- 
tance of ^2200 pc (Diplas & Savage 1994), again with an 
uncertainty of ~ 25%. This places the star on the edge 
of the Sagittarius arm and 250 pc below the plane of the 



Milky Way, well out of the thin disk (scale height <~ 100 
pc, Dickey & Lockman 1990). It is likely that most of 
the gas along this sight line lies within ~ 1000 pc of the 
Sun. This is supported by the fact that the H I column 
toward HD 191877 found by Diplas & Savage (1994) (log 
7V(HI)=20.9) is similar to that toward HD 195965, even 
though the distance to HD 191877 is more than twice as 
great as the distance to HD 195965. The two sight lines 
therefore likely probe a similar distance range in the disk. 
The Hipparcos distance to HD 191877 is d> 1160 pc (ESA 
1997). 

The ground-based optical spectra of HD 195965 and 
HD 191877 in Figure 3 reveal that there are multiple clouds 
along the line of sight to each star. Along both sight 
lines there appear to be two major components of K I, 
as well as several weaker components. The STIS spec- 
trum of HD 195965 (Figure 4) confirms the presence of 
multiple velocity components. Two components are visi- 
ble in the O I 1355.598 A profiles, and the S I 1425.030 A 
profile shows that the strongest of these is actually two 
components. The structure in the S I line matches well 
the K I component structure. The CI I line shows two ad- 
ditional components at more negative velocities, only one 
of which is strong enough to appear in the O I 1355.598 A 
line. The CI I line is usually a good tracer of H 2 (Jura 
1974), so it appears that the molecular hydrogen velocity 
structure along the sight line is not complex. The strong 
O I 1302.168 A line is completely saturated and shows ab- 
sorption extending over ±30 kms -1 , along with damping 
wings at even higher velocities. 

Both sight lines contain a significant amount of H2. 
Most of the unlabeled features in Figures 1 and 2 are H2 
absorption lines. Lines from the lower rotational levels 
(J = 0,1,2) are mostly saturated. These contain most of 
the H 2 column density, so we did not attempt to determine 
the H2 abundance. However, based on the CI I column to- 
ward HD 195965, we estimate roughly (5 — 10) x 10 19 cm~ 2 
on that sight line (Jura 1974). The HD 191877 sight line 
probably contains a similar amount. There are clearly 
multiple components along the HD 195965 sight line (re- 
vealed by the CI I line), but they fall within ±15 kms -1 , 
so at FUSE resolution the H 2 distribution can be treated 
as a single component (see § 4.1). 

The widths of the K I lines in the spectra of both stars 
(and the STIS lines in the HD 195965 spectrum) are larger 
than the resolution of the data. This suggests that there 
are no strong, narrow (unresolved) features on either sight 
line. We can use the STIS data for HD 195965 to show this 
more conclusively for that sight line. Figure 5 shows the 
apparent column density of S I derived from two different 
lines: 1295.653 A and 1296.174 A. Although the oscillator 
strengths differ by a factor of ^3 (Morton 1991), the de- 
rived column density is very similar for the two lines, sup- 
porting the conclusion that narrow, saturated features are 
not present toward HD 195965. (See Savage & Scmbach 
1991 for a discussion of apparent column density profile 
comparisons.) In other words, while cold clouds are clearly 
present, they do not contain a significant amount of the 
material along these sight lines. Most of the column arises 
in the diffuse ISM. 

To test the effects of the lower resolution of FUSE on 
the analysis of blended lines, we smoothed the absorption 
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profiles of the S I 1295.653 A and 1296.174 A lines in the 
STIS data by 15 kms -1 to approximate FUSE resolu- 
tion. We then determined the column densities by sum- 
ming the apparent optical depth over the line profile. The 
total column density derived from the 1295.653 A line is 
N(S I)=(1.21±0.02) x 10 13 cm" 2 , while that derived from 
the 1296.174 A line is N(S I)=(1.34 ± 0.05) x 10 13 cm" 2 . 
Both values are well within the la uncertainties of the 
column densities derived from the unsmoothed profiles, 
N(S I)=(1.26 ± 0.08) x 10 13 cm" 2 from the 1295.653 A 
line, and N(S I)=(1.33 ± 0.18) x 10 13 cm" 2 from the 
1296.174 A line. The measured equivalent widths are 
17.2 ± 0.9 mA and 7.4 ± 1.0 mA, similar to the weaker 
N I and O I lines along the HD 195965 sight line. We 
also carried out this test on the K I lines shown in Figure 
3. Toward HD 191877, the column densities from the un- 
smoothed profile, (4.02 ± 0.06) x 10 11 cm -2 , agrees well 
with the column density derived from the smoothed pro- 
file, (3.97 ± 0.09) x 10 11 cm" 2 . The test was harder to 
perform on the HD 195965 sight line because of residuals 
from poorly removed airglow lines, but the column density 
derived from the unsmoothed profile, (3.24 ± 0.16) x 10 11 
cm.- 2 , is in marginal agreement with that derived from the 
smoothed profile, (2.91 ± 0.19) x 10 11 cm" 2 . With mea- 
sured equivalent widths of 49. 7 ±2.1 mA and 6 2.0 ±1.0 mA 
toward HD 195965 and HD 191877, respectively, these 
lines are similar in strength to the stronger O I, N I, and 
D I lines along both sight lines. These tests provide more 
evidence that there are no unresolved, saturated compo- 
nents along either sight line. 

4. THE D I COLUMN DENSITY 

The simplest approach to analyzing the chemical com- 
position along a line of sight is to measure the equivalent 
widths of all lines of a species, assume the material lies in 
a single component, and construct a Doppler-broadened 
curve of growth (COG). Since the FUSE bandpass con- 
tains all of the transitions of D I except for Lya, in prin- 
ciple this is possible using the current data. In reality, 
many of the D I transitions are not usable due to inter- 
ference from neighboring absorption lines. Table 3 lists 
some of the important interfering lines for the unusable 
D I transitions along these sight lines. The lower order 
D I lines are blended with strong H I lines, while higher or- 
der lines are often blended with strong H 2 lines. However, 
in the spectra of each sight line there are four D I transi- 
tions that were either free of blending or were blended with 
weak H 2 lines that could be removed: A919.102, A920.713, 
A922.900, and A925.974. A more complex but potentially 
more accurate method is to model the absorption lines 
with a profile fitting program. This technique can be es- 
pecially advantageous when the velocity structure of the 
cloud components is complicated. In this section we de- 
scribe both of these approaches to determining the D I 
abundances toward HD 195965 and HD 191877. 

4.1. Molecular Hydrogen Removal 

There are many H 2 lines in the spectral range around 
the deuterium lines, and several can potentially affect the 
measurement of the D I equivalent widths. The lines that 
affect the four D I lines used in this analysis are listed in 
Table 4. We include only lines with rotational level J<5, 



because no higher J level lines were found in the spec- 
tra. The lines listed in Table 4 are expected to be strong 
enough to affect the D I lines. The 920.803 A line does 
not fall within the 920.713 A D I line, but it could affect 
the continuum on the red side of the line. 

To remove the contamination, a line of the same rota- 
tional level in a clean part of the spectrum was found to 
use as a template, following Jenkins et al. (1999). The 
template line was chosen from the same detector segment 
as the line being corrected, in order to avoid differences 
in the line spread function between segments. The pa- 
rameters of the template lines are listed in Table 4. A 
Gaussian function was fitted to the template line, and the 
fit was scaled in optical depth space by the ratio of the /A 
values to have the strength appropriate for the contami- 
nating line. The template line was chosen to have as close 
to the same value of /A as the contaminating line as possi- 
ble in order to avoid saturation effects. The scaled fit was 
then shifted to the wavelength of the contaminating line 
and divided out of the observed spectrum around the D I 
lines. Another H 2 line close in wavelength to the D I line 
was used to register the position of the scaled line before 
dividing it out, in order to compensate for any relative 
errors in the wavelength scale. Figures 6 and 7 show the 
corrected profiles (solid lines) along with the uncorrected 
profiles (dashed lines) for both sight lines as a function of 
LSR velocity. 

4.2. Curve of Growth Analysis 

Once the H 2 lines were removed, it was necessary to 
determine the location of the continuum at the position 
of the D I absorption lines. This was done by fitting low 
order polynomials across the lines. The situation is com- 
plicated by the fact that the D I lines fall on the blue 
damping wings of the associated H I lines. We used the 
steeply-sloped blue sides of the H I lines to constrain the 
fits, and when possible used the red wings of the H I lines 
as a guide in determining the optimal continuum. Figures 
6 and 7 show the absorption profiles of the four D I lines 
after correction for H 2 for HD 195965 and HD 191877, re- 
spectively, along with the adopted continuum placement. 

We measured the equivalent widths of the D I lines us- 
ing the procedures described by Sembach & Savage (1992). 
The lines were integrated over a velocity range that de- 
pended upon the details of the local continuum and line 
blending. The equivalent widths are listed in Table 5 for 
HD 195965 and Table 6 for HD 191877. A theoretical 
single-component COG was fit to the equivalent widths 
by minimizing y 2 while varying the column density N and 
the Doppler parameter b. The final COGs are shown in 
Figures 8 and 9 for HD 195965 and HD 191877, respec- 
tively. The final uncertainties in the equivalent widths 
were determined from the dispersion of the points about 
the COG, assuming that a single component COG is ap- 
propriate for the sight lines. These error bars are a com- 
bination of statistical uncertainties and continuum fitting 
errors, as described by Sembach & Savage (1992), and in- 
clude the effects of fixed-pattern noise. 

Placing the continuum near the 925.974 A line was prob- 
lematic in the spectra of both stars, because of absorption 
from other elements on the blue side of the absorption line 
(see Figures 6 and 7) . This results in large uncertainties in 
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the equivalent widths. Because of this the 925.974 A line 
is not given as much weight as the other lines in the COG 
fitting. Table 7 contains the COG D I column density 
results for both sight lines. 

The inset in each COG figure shows a contour plot of 
the A% 2 values in the N and b dimensions, with contours 
of the minimum A% 2 +1, 4, and 9 shown (corresponding 
to 1, 2, and 3a). The quoted error bars for N shown in 
Table 7 are the maximum extent of the projection of the 
2cr contour onto the N axis, and likewise for b (listed in 
the captions to Figures 8 and 9). 

4.3. Absorption Profile Fitting 

Column densities were also measured independently 
with the profile fitting procedure Owens. f, following the 
method described by Hebrard et al. (2002). Briefly, the 
software simultaneously fits Voigt profiles to all of the un- 
saturated interstellar lines detected in all FUSE segments, 
using a x 2 minimization procedure with many free param- 
eters. The free parameters characterize the physical prop- 
erties of the interstellar clouds (radial velocities, column 
densities, temperatures, micro-turbulent velocities) as well 
as the continuum shapes (modeled with polynomials) and 
possible instrumental effects (line spread function widths, 
wavelength solution variations). The error bars are ob- 
tained using the standard A\ 2 method; they include sta- 
tistical as well as systematic uncertainties (Hebrard et al. 
2002). 

The fits include the D I lines 916.6 A, 919.1 A, and 
920.7 A, the N I lines 951.1 A, 951.3 A, and 959.5 A, as 
well as Fe II, and molecular H2, HD, and CO lines. No 
O I lines were included in the fits because they are all sat- 
urated. Since the H 2 lines were included in the fits, the 
lines were not removed as for the COG analysis (§ 4.1). 
Samples of the fits are plotted on Figure 10, which presents 
two spectral windows for each target. The final fits include 
30 and 22 spectral windows for HD195965 and HD191877, 
respectively. The results for N(D I) and iV(N I) are re- 
ported in Table 7. We obtain logiV(Fe II) = 14.81 ± 0.02 
and 14.95 ± 0.04 (2a) for HD 195965 and HD 191877, re- 
spectively. 

5. THE O I AND N I COLUMN DENSITIES 

We measured the equivalent widths of O I and N I lines 
in the FUSE spectra using the same techniques as de- 
scribed for D I. The measured equivalent widths arc listed 
in Table 5 for HD 195965 and Table 6 for HD 191877. 
Figure 11 shows single-component COGs for O I and N I 
toward HD 195965. The O I lines in the FUSE bandpass 
that are not affected by H2 absorption do not cover a large 
range of oscillator strengths, and they mostly fall on the 
flat part of the COG. The STIS spectrum of HD 195965 
adds two O I lines: the 1355.598 A line on the linear part of 
the COG, and the 1302.169 A line, which exhibits damping 
wings and is on the square root portion of the COG. The 
STIS lines cover a larger range of oscillator strengths, so 
the O I abundance toward HD 195965 is well constrained 
by the FUSE+STIS COG. The N I lines in the FUSE 
spectra also cover a large range of oscillator strengths, so 
an accurate determination of the column density is possi- 
ble. The derived column densities for O I and N I toward 
HD 195965 are listed in Table 7. 



Figure 12 shows the COGs for O I and N I toward 
HD 191877. As was the case for HD 195965, the N I abun- 
dance is well constrained, but the O I abundance toward 
HD 191877 is less certain because of the lack of STIS data. 
A single-component COG fit to the FUSE points yields a 
column density of log N(0 I)= 17.24=bg;|| cm" 2 and the 
Dopplcr broadening parameter b = 8.2+qJ kms _1 (2<7). 
We also set a lower limit of 230 mA (3a) on the 1302.169 A 
equivalent width using an IUE spectrum of HD 191877 
(see § 6). This point rules out the COG found for the 
FUSE points alone, and restricts the O I column density 
to log N(0 I)> 17.42 cm" 2 , with b = 7.85 kms" 1 . The 
6-values derived from the COGs for N I and O I toward 
HD 195965 are almost identical, suggesting that the veloc- 
ity behavior of these two species is coupled. If we assume 
that this is also true for the HD 191877 sight line and 
force b = 7.5 kms" 1 for the O I COG (the solid line in 
the top panel of Figure 12), the derived column density 
of O I is log N(0 I)=17.54±8; 2 £ cm" 2 (2a). Although 
the HD 195965 sight line suggests that O I and N I are 
coupled, we cannot know with certainty whether this is 
also true along the HD 191877 sight line. This exercise 
illustrates the fact that the O I column density toward 
HD 191877 is not well constrained by the lines in the FUSE 
spectrum. We adopt the larger N(0 I) column density as 
our preferred value, but note that it should be treated with 
caution. 

We also used the damping wing s of the O I 1302.169 A 
line in the STIS spectrum of HD 195965 to make an inde- 
pendent determination of the O I abundance toward the 
star. We fit a single component Voigt profile to the nor- 
malized O I line using an oscillator strength f = 0.04887 
and a natural damping constant 7 = 5.75 x 10 8 s" 1 (Mor- 
ton 1991). The best fit model has a central velocity 
Vhelio = +2.1 kms" 1 , b = 6.3±0.3 kms" 1 , andlogA(0 I) 
= 17.80±o;Qg (2a). This column density is in excellent 
agreement with the COG result listed in Table 7. The 
single-component fit to the O I 1302.169 A line is shown 
in Figure 13. The fit shown in the figure reproduces the 
observed absorption profile core and damping wings well, 
with the minor exception of a slight underestimate in the 
absorption strength between +21 and +40 kms -1 . This 
additional absorption is likely due to very low column den- 
sity gas along the sight line that does not contribute much 
to the total column density; a component with an O I 
column density of ~ 10 13 cm -2 could account for this mi- 
nor discrepancy. Experimentation with additional compo- 
nents does not change the best-fit result substantially. Per- 
forming a fit with the O I column density divided among 
two components separated by ~ 12 kms -1 as suggested 
by the O I 1355.598 A line shown in Figure 4 yields a best 
fit nearly indistinguishable from the single-component fit 
shown in Figure 13. 

The O I 1302.169 A line is very strong, and can be used 
to test for the presence of high-velocity components that 
might contribute to the column density of H I but be out- 
side of velocity integration range for other elements. The 
most obvious candidate is the component between +21 and 
+40 kms -1 , which contributes only N(0 I)~ 10 13 cm -2 , 
corresponding to N(R I)~ 3 x 10 16 cm" 2 using the O/H 
ratio of 3.43 x 10" 4 from Meyer (2001), or - 3 x 10" 5 of 
the total H I column density (see § 6). No other compo- 
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nents this strong are visible in Figure 13, so we conclude 
that high-velocity H I is not an issue for the HD 195965 
sight line. 

The STIS spectrum of HD 195965 contains both the op- 
tically thin 1355.598 A line and the damped 1302.169 A 
line. We can use this fact to test for consistency between 
the oscillator strengths of the two lines (see Sofia, Cardclli, 
& Savage 1994). The O I column densities derived from 
these lines are N(0 I)= (5.7 ± 1.1) x 10 17 cm" 2 and 
(6.3±0.5) x 10 17 cm" 2 , respectively (uncertainties are la). 
These results are identical within the uncertainty, which 
suggests that any error the oscillator strength of either line 
is less than 20% (the uncertainty in the 1355.598 A equiv- 
alent width). The oscillator strength for the 1355.598 A 
line listed in Table 5 includes the correction from Welty et 
al. (1999). 

6. THE H I COLUMN DENSITY 

We have re-measured the H I column densities for the 
HD 191877 and HD 195965 sight lines using existing STIS 
and IUE data. Single-component Voigt profiles were fit to 
the H I absorption toward each star. The fits were lim- 
ited to the 1170 to 1260 A spectral region, with the over- 
all continuum levels determined at longer wavelengths, as 
needed. For both sight lines, the fit is insensitive to the 
assumption of a single absorbing component. Adding ad- 
ditional components or changing the 6-value of the ab- 
sorption (assumed to be in the 1 to 100 kms" 1 range) did 
not alter the fits significantly since the quality of the fits 
is governed primarily by the steep walls of the line cores 
and the broad radiation-damping wings. The Voigt profile 
fitting method has been described in detail for H I deter- 
minations in previous FUSE D/H papers (see Moos et al. 
2002, and references therein). Our analyses follow many 
of these same principles and methods. 

For HD 195965, we used the STIS E140H data avail- 
able in this wavelength region (see §2.2) to determine the 
H I column density. We paid special attention to making 
sure that the flux calibration mapping at the ends of the 
individual echelle orders was consistent between adjoining 
orders. Scattered light is not a significant source of uncer- 
tainty in the final measurement since we used the damping 
wings of the lines rather than the saturated portion of the 
line core to determine the H I column density. We find a 
best-fit value of log iV(H I) = 20.95±0.05 (2a) for the fit to 
the high quality STIS data shown in Figure 14. This value 
is slightly higher than, but within the errors of, the for- 
mer value of 20.90+/-0.09 (la) found by Diplas & Savage 
(1994) using the "continuum reconstruction" technique on 
existing IUE data. The reconstruction method and profile 
fitting methods are complimentary to each other. The fit 
is shown in the bottom panel of Figure 14, together with 
the fits appropriate for the ±2cr column density ranges 
(dashed lines). For this sight line, the higher quality STIS 
data more tightly defines the extent of the H I absorption 
than the IUE data, and as can be seen in the top panel 
of Figure 14, the resulting fit is consistent with the IUE 
small aperture data analyzed previously. 

The primary difficulty in fitting the STIS profile of 
HD 195965 arises in the wavelength region between 1222 
and 1225 A, where the fit slightly overestimates the 
amount of absorption present. We explored other mod- 



els that fit this region better (at the expense of fit quality 
at other wavelengths) and found that in all cases the re- 
sulting column density was within the listed uncertainties. 
Therefore, while we are uncertain of the exact cause of 
the discrepancy in this wavelength region, we believe the 
resulting H I column density is a robust estimate. Note 
that this portion of the STIS spectrum contains the bor- 
der between two echelle orders, and that the anomaly is 
not apparent in the IUE spectrum. 

For HD 191877, we re-examined the IUE small aperture 
data (SWP02837) available for this sight line. A more re- 
cent exposure (SWP148425) taken through the small aper- 
ture clearly has problems with the flux levels as a result 
of poor centering; we did not use this low quality spec- 
trum. We find log N(H I) = 21.05±0.10 (2a), which again 
is slightly larger than, but consistent with, the previous 
estimate of 20.90+/-0.10 (la) found by Diplas & Savage 
(1994). The fit is shown in Figure 15, along with the ±2cr 
error bars. As can be seen from this figure, a fit with a 
smaller column, like the one preferred by the previous in- 
vestigators, is clearly below the true column density (see 
the —2a [log N = 20.95] model in the figure). Diplas 
& Savage note that there may be a small amount (<~ 0.01 
dex) of contamination by stellar absorption in the H I Lya 
profile. This contamination does not affect our N(H I) es- 
timate significantly. 

7. DISCUSSION 

The sight lines analyzed in this paper are the longest 
ones for which the D I column densities have been deter- 
mined with FUSE. As such, they present the opportunity 
to look for variations in D I/H I, D I/O I, and D I/N I 
beyond the local ISM sight lines presented by Moos et al. 
(2002). Table 8 contains column density ratios of various 
elements toward both stars, along with the mean values for 
the local ISM from Moos et al. (2002). In this section we 
discuss the behavior of these column density ratios beyond 
the local ISM. 

7.1. D/H 

The deuterium abundances along seven sight lines have 
been analyzed with FUSE (Moos et al. 2002; Kruk et al. 
2002; Friedman et al. 2002; Sonneborn et al. 2002; Lemoine 
et al. 2002; Lehner et al. 2002; Wood et al. 2002; Hebrard 
et al. 2002). Five of these had both D I and H I mea- 
surements, and the results were consistent with a constant 
value of D /H in the local ISM. Earlier Copernicus obser- 
vations showed evidence for variations in D/H at larger 
distances (see Vidal-Madjar & Gry 1984 for a review), ev- 
idence which was strengthened by subsequent IMAPS in- 
vestigations (Jenkins et al. 1999; Sonneborn et al. 2000). 
The main obstacle to verifying and understanding this 
variation is the paucity of suitable deuterium sight lines 
that probe distances >250 pc. 

Figure 16 shows the D/H ratios for HD 195965 and 
HD 191877 plotted as a function of distance from the Sun, 
along with five of the seven original FUSE sight lines with 
reliable H I measurements, and several others from HST, 
IMAPS, and Copernicus (sec Moos et al. 2002). The D/H 
ratios toward both stars are lower than the average value 
found for the original seven FUSE sight lines. This con- 
trasts with the lack of variation seen in the local ISM. The 
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new sight lines support the existence of variations in the 
D/H ratio at distances larger than a few hundred pc, as 
also suggested by IMAPS and Copernicus. 

In the simplest model of chemical evolution, stellar pro- 
cessing (astration) leads to a decline in the D/H ratio be- 
cause D is easily destroyed in stars. Thus, the D/H ratio 
is expected to vary among regions of the Galaxy that have 
undergone different amounts of star formation in the past, 
provided that the gas in these regions has not mixed {e.g. 
Tosi et al. 1998; Chiappini et al. 2002). In a small volume 
like the local ISM (r = 100 pc), Moos et al. (2002) showed 
that the time between supernova tsN is roughly equal to 
the mixing time scale (assumed to be the sound crossing 
time t s ), suggesting that the gas in such a region is well- 
mixed. If a larger region is considered, say 1000 pc in ra- 
dius to approximate the scales probed by HD 195965 and 
HD 191877, the crossing time increases by a factor of 10, 
while the time between supernovae decreases by a factor 
of 100, so that tsN -C t s . Supernovae are not distributed 
randomly over this region but are spatially correlated in re- 
gions of star formation, meaning that substantial inhomo- 
geneities will arise between star-forming regions and more 
quiescent regions. For these reasons it is not surprising 
to see variations over the distances probed by HD 195965 
and HD 191877. 

A more precise treatment of mixing in the ISM was car- 
ried out by de Avillez & Mac Low (2002). Using numerical 
simulations, they found that the mixing time scales were 
a few times 10 8 years. Chemical inhomogeneities in the 
ISM are long-lived, so variations in D /H are expected over 
1000 pc length scales. In fact, since the mixing time scales 
increase only slightly with length scale in their models, it is 
unlikely that the local ISM (e.g., the original seven FUSE 
sight lines) would have a constant D/H value if any sources 
of inhomogeneities (i.e., supernovae) have occurred within 
the past ~ 10 8 years. 

We have investigated several alternative explanations 
for the low D/H values. One possibility is that some of 
the D column is in high-velocity components which are 
outside of the velocity integration range used to deter- 
mine N(D I) (Jenkins et al. 1999). We have used the O I 
1302.169 A line in the STIS spectrum of HD 195965 to 
rule this out for that sight line (see § 5). Note that <~ 50% 
of the D I column would have to be missed if the D /H ra- 
tio for the HD 191877 sight line were as high as the Moos 
et al. (2002) mean local ISM value. Figure 7 shows that 
this is not the case, although we cannot rule out a much 
smaller systematic offset for HD 191877. 

A second possibility is that a significant amount of D 
is locked up in HD molecules. To test this idea, we 
have measured the equivalent widths of several strong, 
isolated Lyman series HD lines in the FUSE spectra: 
(3-0) R(0) 1066.271 A, (5-0) R(0) 1042.847 A, and (4- 
0) R(0) 1054.288 A. The column densities derived from the 
equivalent widths are ~ 2.4 x 10 14 cm -2 for the HD 195965 
sight line, and - 3.0 x 10 14 cm~ 2 for the HD 191877 sight 
line. For both sight lines this is < 3.5% of the total D col- 
umn density, much lower than the uncertainty in N(T) I). 

7.2. O/H 

The D/O and O/H ratios toward HD 195965 in Table 8 
are quite different from the average values in the local ISM 



found by Moos et al. (2002). This is a surprising result, 
since Meyer et al. (1998) found very little variation in O/H 
in their survey, which covered distances and H I column 
densities similar to those of the sight lines in this paper. 
Both the O I and H I column densities toward HD 195965 
were determined through at least two independent meth- 
ods, with very similar results. The H I was measured using 
both IUE and STIS data, and the O I was measured with a 
single-component COG which combined FUSE and STIS 
data, as well as independent estimates from the optically 
thin 1355.598 A line and the damped 1302.169 A line in 
the STIS spectrum. Both of these column densities are 
therefore quite firmly established, and we believe that it is 
very unlikely that there is a substantial systematic error 
in either quantity. 

HD 195965 appears to lie on a very high metallicity sight 
line, with supcrsolar O/H (Sofia & Meyer 2001). This 
may be a region where the signature of stellar processed 
gas is still visible (as is not the case for the local bubble). 
Humphreys (1978) found that this star is in the Cygnus 
OB7 association. If so, it may be in a peculiar region, per- 
haps in a metal-enriched cloud of ejecta from the star itself 
or a neighbor. Since gas with high metallicity is thought to 
have undergone stellar processing, D/H in such gas is ex- 
pected to be low. This expected anticorrelation between 
O/H (or N/H) and D/H was not seen in the first seven 
FUSE sight lines, although it could have been masked by 
the small number of sight lines studied (Stcigman 2002). 
Figure 17 compares the O/H and D/H value for HD 195965 
with those of the original seven FUSE sight lines. Unfor- 
tunately HD 195965 is the only point at high metallicity, 
so no firm conclusions can be drawn. However, it is pos- 
sible that this sight line shows evidence of the expected 
anticorrelation in the ISM (see also Steigman 2002). Note 
that an error in the H I column density toward HD 195965 
would tend to move that point along a diagonal line with 
a positive slope, and thus could not produce the expected 
trend of decreasing D/H with rising O/H. 

The HD 191877 sight line lacks STIS data, so we could 
not constrain the O I and H I abundances with multiple 
estimates. For these reasons we are less confident in the 
O I column densities found for the HD 191877 sight line. 
As discussed in § 5, the derived log N(0 I) is very un- 
certain, although our preferred value, found by coupling 
the 6-values of N I and O I, is consistent with the FUSE 
local ISM ratio and the Meyer et al. (1998) ratio. Unfor- 
tunately, with the existing information, we cannot make 
any strong conclusions about O/H toward HD 191877. 

7.3. N/H 

The N/H ratios toward HD 195965 and HD 191877 in 
Table 8 are higher than the average values in the local 
ISM found by Moos et al. (2002). Figure 18 shows a plot 
of N/H versus D/H along both sight lines, compared to 
the four FUSE sight lines for which N/H was measured. 
The figure is very suggestive of an anticorrelation between 
D/H and N/H. ' 

Complicating the interpretation of this plot is the fact 
N/H ratio is susceptible to ionization effects, especially 
at low column densities which do not provide adequate 
shielded from ionizing photons, such as the column densi- 
ties along the local ISM sight lines (Jenkins et al. 2000). 
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The dashed line in Figure 18 shows the ISM value of N/H 
found by Meyer et al. (1997). The sight lines used in that 
study had higher H I column densities than the FUSE lo- 
cal ISM sight lines, so the FUSE sight lines are very likely 
underestimating N/H, because a significant fraction of N 
is ionized. The apparent trend in Figure 18 may simply 
reflect the tendency for N(N II)/iV(N I) to decrease with 
increasing ]V(H I). 

However, an anticorrelation between D and N abun- 
dances would be expected if D is destroyed in stars while 
N is produced (e.g. Steigman et al. 2002). The over- 
abundance of heavy elements accompanied by the under- 
abundance of D seems to fit well within a simple chemical 
evolution scenario. While ionization effects probably play 
a role in the high N/H ratio on these sight lines, they may 
not be the only factor. A comprehensive N/H survey using 
FUSE spectra is currently underway, and may shed light 
on this issue. 

Several sight lines do not follow the expected trend of 
low D/H accompanied by high N/H or O/H. Sonneborn 
et al. (2002) found a low O/H ratio toward BD+28°4211 
with FUSE, while their D/H ratio is close to the Local 
ISM mean. Jenkins et al. (2000) found a low D/H toward 
5 Orionis using IMAPS, but also a somewhat low N/H 
ratio. These deviations suggest that the simple picture 
of increasing heavy elements abundances with decreasing 
D/H is not complete. 

7.4. Future Prospects for D/H Measurements 

Measuring D/H values in the Milky Way becomes in- 
creasingly difficult as N(R I) increases, perhaps reaching 
a fundamental confusion limit. The COGs in Figures 8 
and 9 show that the high order D I Lyman transitions 
are close to saturation. At higher column densities, these 
transitions will fall on the flat part of the COG and will 
not give precise column density information. Even higher 
order transitions at shorter wavelengths, where the D I 
lines might not yet be saturated at higher column den- 
sity, are not useful because the H I transitions begin to 
blend together. These facts make it unlikely that many 



Galactic disk sight lines more distant than those studied 
here will be useful for D/H measurements. However, the 
sight lines analyzed in this paper demonstrate that D/H 
properties over the distance range from 500 to 2500 pc can 
be explored with FUSE, even if these sight lines contain 
a significant amount of H 2 . This first glimpse at D/H in 
this distance range suggests that it will be very useful for 
understanding mixing and chemical evolution in the ISM. 

The chemical evolution model of Chiappini et al. (2002) 
predicts an increase in D/H with galactocentric radius, 
accompanied by a decrease in N/H and O/H, leading to 
a pronounced gradient in D/O and N/O. The gas in the 
outer parts of the galaxy has undergone less stellar pro- 
cessing, and in fact Chiappini et al. (2002) predict that the 
D/H ratios in the outer disk will approach the primordial 
value. The sight lines studied here have roughly the same 
galactocentric radii as the local ISM sight lines of Moos et 
al. (2002), so we cannot test this model, but with FUSE 
observations of sight lines in the inner and outer Galactic 
disk this may prove possible. 
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Fig. 1. — A portion of the FUSE SiC2 spectrum of HD 195965. The Lyman transitions of H I and D I are labeled, as are the O I and N 
I transitions in this region of the spectrum that were used in the analysis. The D I transitions used in the curves of growth were Ly7, Ly8, 
Ly9, and LylO. Most of the other features are H2 absorption lines. The spectrum shown here is binned by 5 pixels to ~10 km s _1 bins for 
display. 
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Fig. 2. — A portion of the FUSE SiC2 spectrum of HD 191877. The Lyman transitions of H I and D I are labeled, as are the O I and N 
I transitions in this region of the spectrum that were used in the analysis. The D I transitions used in the curves of growth were Ly7, Ly8, 
Ly9, and LylO. Most of the other features are H2 absorption lines. The spectrum shown here is binned by 5 pixels to ~10 km s _1 bins for 
display. 
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Fig. 3. — Absorption profiles of the K I 7698.964 A line for both stars. The absorption lines at -20 kms and +28 kms are remnants 
of atmospheric H2O lines. 
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Fig. 4. — Absorption profiles of selected lines from the STIS spectrum of HD 195965. 
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Fig. 5. — Comparison of the apparent column density profiles derived from two S I lines in the STIS spectrum of HD 195965. The strengths 
(/A) of the two lines differ by a factor ~3, but the apparent column densities are similar at all velocities. The total column density between 
3 and 13 km s _1 derived from the 1295.653 A line is N(S I)=(1.26 ± 0.08) X 10 13 cm" 2 , while that derived from the 1296.174 A line is 
N(S I)=(1.33 ± 0.18) X 10 13 cm -2 . A representative error bar is shown. 
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Fig. 6.— Four D I absorption profiles in the FUSE HD 195965 spectra. The 919.102 A, 920.713 A, and 922.900 A profiles have been 
corrected for H2 absorption. The uncorrected profiles are also shown (dashed lines). The 925.974 A profile is not affected by H2 absorption. 
The adopted continua are shown as dotted lines. The feature at 925.780 A (—45 kms -1 in the bottom panels) is unidentified. 
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Fig. 7.— Four D I absorption profiles in the FUSE HD 191877 spectra. The 919.102 A, 920.713 A, and 922.900 A profiles have been 
corrected for H2 absorption. The uncorrected profiles are also shown (dashed lines). The 925.974 A profile is not affected by H2 absorption. 
The adopted continua are shown as dotted lines. The feature at 925.780 A (—45 kms -1 in the bottom panels) is unidentified. 
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Fig. 8.— The best fit single-component curve of growth for D I toward HD 195965, which gives log jV(D I)=15.83±g;}f (2a) and b = U.1±IY 
km s _1 (2a). The filled circles are points measured from the SiCl channel, and the open circles are points measured from the SiC2 channel. 
The inset shows a contour plot of the \ 2 distribution, with the la, 2a, and 3a contours shown. 
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Fig. 9.— The best fit single-component curve of growth for D I toward HD 191877, which gives log iV(D I)=15.95±g;f£ (2a) and b = 9.9±4.0 
km s — 1 (2<j). The filled circles are points measured from the SiCl channel, and the open circles are points measured from the SiC2 channel. 
The inset shows a contour plot of the \ 2 distribution, with the lu, 2<r, and 3<r contours shown. 
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Fig. 10. — Examples of Owens. f absorption profile fits. The histogram lines are the data, the solid thick lines are the fits, and the dotted 
lines are the continua. Absorption features are labeled. 
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Fig. 11. — The best fit single-component curve of growth for O I (top panel) and N I (bottom panel) toward HD 195965. The COG 
parameters for O I arc log N(0 I)=17.74±g;gjj (2a) and b = 5.9 ± 0.8 km s" 1 (2a), and the parameters for N I arc log JV(N I)=16.92±g;g| 
(2a) and b = 5.8 ± 0.3 km s" 1 (2a). The filled circles are points measured from either the FUSE SiCl or LiFl channel or from STIS, and the 
open circles are points measured from either the FUSE SiC2 or LiF2 channel. The insets show contour plots of the A\ 2 distributions, with 
the la, 2a, and 3a contours shown. 
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Fig. 12. — The best fit single-component curve of growth for O I (top panel) and N I (bottom panel) toward HD 191877. The COG 
parameters for O I are log JV(0 I)=17.54±g;f§ cm" 2 (2a) if b is forced to be 7.5 km s" 1 as suggested by the N I COG (the COG for these 
parameters is shown as the solid line in the top panel). If b is allowed to vary, we find log N(0 I)=17.24±fj || cm~ 2 (2a) and b = 8.2±q'8| km 
s _1 (2a), the curve for which is shown as a dashed line in the top panel. The lower limit on the 1302.169 A line (the arrow in the top panel) 
restricts the O I column density to log N(0 I)> 17.42 cm . The parameters for N I are log N(N I)=16.85±g;p| (2a) and b = 7.5 ± 0.5 km 
s _1 (2a). The filled circles are points measured from either the FUSE SiCl or LiFl channel, and the open circles are points measured from 
cither the FUSE SiC2 or LiF2 channel. The insets show contour plots of the Ax 2 distributions, with the la, 2a, and 3a contours shown. 
The plus sign in the inset for O I shows where b and N for the solid (b = 7.5 km s— 1) COG fall on the Ax 2 contours. The large error bars 
for some O I equivalent widths are a result of continuum placement uncertainty. 
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Fig. 13. — The continuum-normalized O I 1302.169 A line profile from the STIS spectrum of HD 195965. A single-component model fit to 
the damping wings is also shown. The best fit model gives log N(0 I)=17.80 and b=6.3 kms -1 . 
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Fig. 14. — H I Lya absorption toward HD 195965. The top panel shows small-aperture IUE date from Diplas & Savage (1994). The bottom 
panel show the STIS spectrum. The best fit Voigt profile (solid line) and the ±2<r profiles (dashed lines) arc shown. The model shown is for 
log JV(H I)= 20.95 ± 0.05. 
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Fig. 15. — H I Ly« absorption toward HD 191877. The spectrum is small-aperture IUE date from Diplas & Savage (1994). The best fit 
Voigt profile (solid line) and the ±2<x profiles (dashed lines) are shown. The model shown is for log N(R T)= 21.05 ± 0.10. 
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Fig. 16. — D/H measurements for 17 sight lines, plotted versus distance to the continuum source. The sight lines studied in this paper are 
labeled, as are the sight lines observed by IMAPS and the Capclla sight line. Following Moos et al. (2002), the error bars shown are 2<r/2. 
Adapted from Moos et al. (2002). 
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Fig. 17. — O I/H I for HD 195965, along with five of the seven original FUSE sight lines with reliable H I measurements (Moos ct al. 2002). 
Following Moos ct al. (2002), the error bars shown are 2<r/2. The dashed line marks the O/H value found by Meyer et al. (1998) for the local 
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Table 1 
Sight Line Information 21 



Name 


HD 195965 


HD 191877 


I 


85?71 


61°57 


b 


5°00 


-6°45 


Distance b (pc) 


794 ± 200 


2200 ± 550 


z dist. (pc) 


69 


-249 


V (mag) 


6.98 


6.26 


E(B-V) 


0.25 


0.18 


Spectral Type 


BOV 


Bllb 


Log [N(HI) (cm- 2 )] (IUE) 


20.90 ± 0.09 


20.90 ± 0.10 



a From Diplas & Savage (1994) and references therein. 

b We have assumed 25% uncertainty associated with the 
spectroscopic parallax measurements. 
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Table 2 



Observation Log 



Star 


Instrument /Grating 


Datasct ID 


Observation Date 


Exposure Time 
(s) 


Aperture 


A coverage 

(A) 


HD 195965 


FUSE 


P1028802, P1028803 


2000 Nov 8,9 


30000 


LWRS 


905 - 1187 


HD 195965 


FUSE 


P1028801 a 


2000 Jun 20 


6400 


LWRS 


905 - 1187 


HD 195965 


STIS/E140H 


O6BG01010 


2001 Oct 9 


415 


0.1" x 0.03" 


1197 - 1401 


HD 195965 


IUE 


SWP10843 


1980 Dec 20 


1560 


Small 


1150 - 1970 


HD 191877 


FUSE 


P2051101 


2001 Jul 30 


28100 


LWRS 


905 - 1187 


HD 191877 


FUSE 


P1028701" 


2000 Jun 5 


6100 


LWRS 


905 - 1187 


HD 191877 


IUE 


SWP02837 


1978 Oct 2 


2280 


Small 


1150 - 1970 



a Datascts P1028701 and P1028801 were short snapshots. These data were included for the profile fitting analysis, but not for the 
curves of growth. 





Table 3 




Interfering Lines 


A 


Contaminating Lines 


(A) 




D I 972.272 


H I Ly7, O I A972.143 


D I 949.484 


H I Ly<5, H 2 Lyman (14-0) P(2) A949.353, H 2 Werner (3-0) P(2) A949.611 


D I 937.548 


H 2 Werner (4-0) P(4) A937.557, Fe II A937.652 


D I 930.495 


I A930.257, H 2 Werner (4-0) R(2) A930.446, H 2 Werner (4-0) Q(l) A930.578 



Table 4 

Interfering H 2 Line Removal Parameters 



Interfering H 2 Template H 2 



D I A 


ID 


A 


log A/ a 


ID 




A 


log A/ a 


(A) 




(A) 






(A) 


919.102 


Werner (5 - 0) R(4) 


919.047 


0.931 


Lyman (13 — 


0) R(4) 


962.150 


0.929 


920.713 


Werner (5 - 0) 0(4) 


920.830 


0.0678 


Lyman (12 — 


0) 7?(4) 


970.830 


0.194 


922.900 


Lyman (18 - 0) P(3) 


922.893 


-0.820 


Lyman (17 — 


0) R(4) 


928.440 


-0.511 


925.974 


none 















a Oscillator strengths from Abgrall, Roueff, & Drira 2000. 



Hoopes et al. 



Table 5 

Measured Equivalent Widths Toward HD 195965 



Wavelength 


lncr A f a 


W^ fSiOl or LiFI *) h 


W, fSiO? or FiF?"! 15 

VV \ \ kJlV^Zj Ul J_j1_L Zj ) 




( A ) 








DI 


919.102 


0.0425 


48.8 ± 3.9 


38.4 ± 4.3 


DI 


920.713 


0.171 


54.9 ±5.3 


64.3 ±8.4 


DI 


922.900 


0.312 


66.1 ±12.2 


74.2 ±3.5 


DI 


925.974 


0.469 


77.4 ±24.9 


76.5 ± 10.2 


01 


921.860 


0.0402 


84.9 ±3.6 


90.3 ±3.4 


01 


922.220 


-0.646 


75.1 ±4.2 


70.0 ±2.7 


01 


936.630 


0.527 


93.3 ±3.6 


99.6 ±3.8 


01 


1039.230 


0.980 


120.4 ±4.6 


126.1 ±4.8 


01 


1039.230 


0.980 


115.2 ±4.4 


123.5 ±4.8 


01 


1302.169 c 


1.804 


231.2 ±5.9 




01 


1355.598 c 


-2.803 


10.7 ±2.1 




NI 


951.079 


-0.911 


40.5 ±3.4 


43.9 ±4.1 


NI 


951.295 


-1.766 


12.4 ±2.4 


12.8 ±2.9 


NI 


953.415 


1.100 


92.0 ±4.4 


92.6 ±4.6 


NI 


953.655 


1.377 


93.6 ±4.6 


95.6 ±4.6 


NI 


959.494 


-1.321 


21.3 ± 1.9 


19.6 ±2.2 


NI 


1134.165 d 


1.237 


116.9 ±5.6 




NI 


1134.165 d 


1.237 


113.9 ±5.6 




NI 


1159.817 d 


-2.006 


9.5 ± 1.2 




NI 


1160.937 d 


-2.553 


1.9 ± 1.0 





a Oscillator strengths from Morton 1991, except that for the OI 
1355.598 A line, which is from Welty 1999. 

b Lines with A < 1050 A were measured from the SiC spectra, and 
lines with A > 1000 A were measured from the LiF spectra. 

c These measurements are from the STIS spectrum. 

d Detector x-walk effects due to the bright airglow lines at these 
wavelengths prevented the measurement of the line strengths in the 
LiF2 channel. 
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Table 6 

Measured Equivalent Widths Toward HD 191877 



Wavelength 


log x.r 


W A (SiCl or LiFl) b 


W A (SiC2 or LiF2) b 








(mA) 


' (mA) 


DI 


919.102 


0.0425 


COO f A A 

52.8 ± 4.9 


ATA 1 A CT 

47.4 ± 4.5 


DI 


920.713 


0.171 


G A A 1 11 O 

64. U ± 11.8 


OA A 1 A O 

61.4 ± 4.8 


DI 


922.900 


A O 1 O 

0.312 


no C 1 G A 

08.6 ± 6.4 


£?7 O 1 C O 

67.8 ± 6.8 


DI 


925.974 


A ARf\ 


TO O _L OO 

78. o ± 22.5 


OA 7 I t; /I 

8U.7 ± 5.4 


OI 


921.860 


A A A AO 

(J.U4U2 


AO A 1 A A 

98. U ± 4.4 


1 AO E 1 A 1 

l(Jo.5 ± 4.1 


OI 


922.220 


A C A f 

-0.646 


77.9 ± 6.2 


OA A I o a 

80.0 ± 8.6 


OI 


925.446 


A A AA 

-0.490 


O O E 1 1 A O 

82.5 ± 10.2 


OCO 1 1 1 A 

85.3 ±11.9 


OI 


936.630 


0.527 


119.6 ± 2.4 


1 1 O 1 I O A 

118.1 ± 2.4 


OI 


1039.230 


A AO A 

0.980 


1 A C A I on 

145.0 ± 2.9 


145.1 ± 2.9 


OI 


1039.230 


0.980 


142.6 ±2.9 


142.6 ±2.9 


NI 


951.079 


-0.911 


39.4 ±7.0 


36.9 ±6.7 


NI 


951.295 


-1.766 


10.8 ±4.6 


12.5 ±2.6 


NI 


953.415 


1.100 


118.1 ±6.2 


118.1 ±6.2 


NI 


953.655 


1.377 


120.6 ±6.2 


119.8 ±6.2 


NI 


959.494 


-2.006 


18.0 ±2.8 


20.9 ±3.4 


NI 


1159.817 c 


-2.553 


9.3 ± 1.8 




NI 


1134.165 c 


1.237 


132.2 ±6.7 




NI 


1134.415 c 


1.528 


141.3 ±7.2 




NI 


1134.980 c 


1.693 


159.8 ±8.3 





a Oscillator strengths from Morton 1991. 

b Lines with A < 1050 A were measured from the SiC spectra, and 
lines with A > 1000 A were measured from the LiF spectra. The O I 
1039 A lines was measured in both LiF and SiC channels. 

c Detector x-walk effects due to the bright airglow lines at these 
wavelengths prevented the measurement of the line strengths in the 
LiF2 channel. 



Table 7 
Column Densities' 1 



Name 


COG 




HD 195965 
Profile fit 


Adopted Value b 


HD 

COG 


191877 

Profile fit 


Adopted Value b 


log N(D 1) 
log N(H I)° 
log N(0 l) d 
log N(N I) 


15.83±g 

17.74±g 
16.92±j| 


14 
11 

05 

81 

06 


15.97±°;;| 
20.95 ± 0.05 
17.80±°;°| 
16.83±°;[» 


15.88iQ ^3 
20.95 ± 0.05 
17.77±°|| 
16.85± ; 06 


i5.95±g:?i 

17.54±g 12 
16.85± ; 08 


15.93±°;;? 
21.05 ± 0.10 

16.93±°;?^ 


15.94±°;? 2 
21.05 ± 0.10 
17.54±«?r 
16.88±°;?« 



a Error bars arc la estimates (95% confidence). For all of the adopted values except O I the uncertainties 
were chosen to be the largest uncertainties of the individual values, because the two measurements were made 
on the same data and arc therefore not independent. For O I toward HD 195965 the measurements are largely 
independent, so the errors were added in quadrature. 

b For D I and N I, the adopted value is the weighted mean of the single-component curve of growth and profile 
fitting results. When only one method was available (e.g., H I), the result from that method is the adopted 
value. 

c The H I column densities were derived by model fits to the damping wings of Lya. For HD 195965, Lya was 
observed with STIS, and for HD 191877 we reanalyzed IV E data from Diplas & Savage 1994. 

d The O I column density toward HD 195965 was derived by a single-component curve of growth and a model 
fit to the damping wings of the 1302.169 A line in the STIS data. 

°The O I column density listed in the table is the value derived when b is set to 7.5 km s — 1 . If b is allowed 
to vary, the best fit COG gives log A^(0 I) — 17.24±q;22 cm -2 (see discussion in § 5). 
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Table 8 
Column Density Ratios' 1 



Name HD 195965 HD 191877 Local ISM b 



D/H 


(0.85±g-|t) x 10" 


-5 


(o.78±8:§§) x io- 


-5 


(1.52 ±0.15) x 10" 


-5 


D/O 


(1.29±g-^) x 10" 


-2 


(2.51±2;i4) x 10- 


-2c 


(3.99 ±0.38) x 10" 


-2 


D/N 


(i.07±g:M) x io- 


-1 


(1.15±°; 3 7 I) x 10- 


-1 


(3.30 ±0.40) x 10" 


-1 


0/H 


(6.61±J;-?f) x 10" 


-4 


(3.09±J-||) x 10" 


-4c 


(3.03 ±0.42) x 10" 


-4 


N/H 


(7.94±i-g|) x 10" 


-5 


(6.76±f P) x 10" 


-5 


(4.24 ±0.62) x 10- 


-5 


0/N 


8.32±i;66 




4.57±?;i 




8.1 ±0.8 





a Error bars are 2a estimates (95% confidence). 

b Mean Local ISM values from Moos et al. 2002 and references therein. 
The uncertainties for these values are 2a in the mean. 

c Thc D/O and O/H ratios toward HD 191877 listed in the table use the 
O I column density found by setting b to 7.5 km s _1 . If b is allowed to vary, 
wc find D/O=(5.01±l If) x 10~ 2 , 0/H=(1.55±g-|3) x 10~ 4 (sec discussion 
in § 5). 



